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Abstract 

Using the level-spacing distribution and the total probability function of the 
numbers of levels in a given energy interval we analyze the crossover of the 
level statistics between the delocalized and the localized regimes. By numer- 
ically calculating the electron spectra of systems of up to 32 3 lattice sites de- 
scribed by the Anderson Hamiltonian it is shown that the distribution P(s) of 
neighboring spacings is scale- independent at the metal-insulator transition. 
For large spacings it has a Poisson-like asymptotic form P(s) oc exp(— As/ A), 
where A ~ 1.9. At the critical point we obtain a linear relationship between 
the variance of the number of levels {[5n(e)] 2 ) and their average number (n(e)) 
within the interval e. The constant of proportionality is less than unity due 
to the repulsion of the levels. Both P(s) and ([<5n(e)] 2 ) are determined by the 
probability density Q n (e) of having exactly n levels in the energy interval e. 
The distribution Q n (e) at the critical point is found to be size-independent 
and to obey a Gaussian law near its maximum, where n ~ (n). 



*to be appeared in Jpn. J. Appl. Phys., 34 (8A), 1995. 
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I. INTRODUCTION 



The statistical properties of the energy levels of disordered systems near the metal- 
insulator transition are presently attracting much attention. Due to Anderson localization, 
the one-electron states of systems with higher dimensionality (d > 2) are known to experience 
a crossover from delocalized to localized behavior. This metal-insulator transition (MIT) is 
induced by increasing the fluctuations of a random potential or by shifting the Fermi energy 
through the mobility edge at fixed disorder. M The statistics of the corresponding energy 
levels is very sensitive to the MIT and depends dramatically on the strength of disorder. 
In the insulating regime the levels are completely random, and their distribution obeys the 
Poissonian statistical laws of uncorrelated variables. On the other side of the MIT, in the 
metallic regime, correlations in the spectrum become stronger due to quantum-mechanical 
level repulsion. It was shown earlier |2|,§] that the level statistics in a disordered metal is 
governed by the random-matrix theory. [[§-|6|] 

By using finite-size scaling arguments Shklovskii et al. |7j suggested that exactly at the 
MIT the distribution of spacings P(s) between the neighboring energy levels exhibits critical 
behavior and is a scale-invariant function which differs considerably from both the Wigner 
surmise Pw(s) = (its/ 2) exp(— 7rs 2 /4) and the Poisson law Pp(s) = exp(— s), corresponding 
to the delocalized and the localized regime, respectively. Here s is measured in units of 
the mean level spacing A. Later the universal properties of P(s) at the mobility edge were 
studied in some detail by using approximate analytical ||[§ and numerical methods. [|T0| -[l2 



In this paper we present results of elaborate numerical investigations of the statistics 
of fluctuations in the energy spectra of three-dimensional (3D) disordered lattices near the 
MIT. Using the Anderson model we calculate very accurately not only the critical distri- 
bution P(s) of neighboring spacings, but also the total probability density Q n {e) of having 
n levels within an energy interval e. The latter contains more complete information about 
spectral correlations of higher orders. Both distributions exhibit a universal character at 
the critical disorder and become independent of the system size. We find that the variance 
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of the number of the levels in an interval e as a function of mean number of levels is also 
scale-invariant, and obeys a linear relation (5n 2 ) = n(n) with k « 0.27, when (n) ^> 1. 

II. NEAREST-NEIGHBOR LEVEL SPACING DISTRIBUTION 

In order to calculate the spectrum we use the Hamiltonian corresponding to a tight- 
binding model with diagonal disorder, H = J2 n £ n a n a n + J2 n ^m( a n a m + c.c), where (a n ) 
is the creation (annihilation) operator of an electron at a site n in a simple cubic lattice. The 
random site energy e n is distributed uniformly between — W/2 and W/2. The second sum 
involves only pairs of nearest sites {n,m}. The MIT in the center of the band corresponds 
to W c w 16.5. |3 

The energy levels of cubes of various linear sizes from L = 5 to 32 with periodic boundary 
conditions were calculated by direct diagonalization of the Hamiltonian using an optimized 



Lanczos algorithm. []12| Figure [I] displays the histogram of P(s) calculated at the disorder Fig. 1 



corresponding to the critical point, W c . The number of realizations was chosen so that the 
numbers of levels from central half part of the band were approximately 1.2 10 7 , 1.2 10 6 and 
2 10 5 for L = 5, 12 and 32, respectively. All of the obtained data, regardless of L, lie on the 
same curve, the critical P(s), which is considerably different from both P\v{s) and Pp(s). 
However, for disorder W ^ W c the distribution P(s) depends on the system size and exhibits 
scaling behavior near the MIT. [|7| Indeed, for decreasing W the distribution of spacings 
approaches the Wigner surmise, while for increasing W it scales towards a Poissonian law. 
This crossover is apparently accompanied by the transition from the delocalized to the 
localized regime. The finite-size scaling properties of P(s) allow determination of the critical 
exponent v ~ 1.45 and the disorder dependence of the scaling parameter, the correlation 
length £(W). [0,0 We have also verified that the form of P(s) exactly at the MIT does 



not depend on the width of the energy interval as long as the levels in this interval belong to 
the critical region, i.e. satisfy the condition L < £ = (|e — e c \/e c )~ u , where e c is the mobility 
edge. 
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Using a large number of realizations of the randomness we were able to analyze the 
asymptotic behavior of the critical P(s) as shown in the inset of Fig. |1|. For spacings s > 2 
the exponential tail of P(s) is intermediate between the Gaussian and the Poissonian. It is 
well described by the relation P(s) = exp(— As) with a numerical coefficient A « 1.9. One 
observes that this asymptotic decay deviates from the power-law In P(s) oc s 1+1 ' 3v which 
was recently found || by using the effective "plasma model" as defined by Dyson. J5| 

III. PROBABILITY OF N LEVELS IN A GIVEN INTERVAL 

In order to analyze the correlations between several consecutive eigenvalues, one can 
study the probability that an energy interval of width e centered at a randomly chosen 
energy contains exactly n levels. The distribution of this joint probability Q n (e) is related 
to the n-level correlation function R n (ei,E2, ...,e n ). It provides a more complete description 
of the level statistics than the two-point correlation function, i?2(£i — £2), which can 
be expressed in terms of the n-level spacing distribution p(n, s), R^s) = J2^=oP( n , s )- On 
the other hand, the probability of finding no levels (n = 0) inside the interval defines the 
nearest-neighbor spacing distribution, P(s) = d 2 Qo(s)/ds 2 . 

The results of the random-matrix theory for Q n (e), which correspond to the metallic 
case, have been calculated numerically, || and were later expressed in an explicit analytical 
form |l4j for the three universality classes of random Hamiltonians: orthogonal, unitary and 
symplectic. In the limit 5n = \n — e\ <C s the distribution for the orthogonal ensemble is 



approximately Gaussian, 



where e is measured in A and B depends weakly on Sn. 

In the insulating regime the sequence of the levels is completely random due to the 
localization of the states, and therefore we have the usual Poissonian process, 

Q n (e) =e n exp(-e)/n\. (2) 



For example, if n = the Wigner surmise and the Poisson law give 



Qo(e) = l~4= Err (^)> Qo(e) = exp(-e) (3) 

'71 I 



for the metallic and insulating regime, respectively. Err(x) is the error function. 

It is of great interest to study how Q n (e) changes from ([!]) to (0) with increasing disorder 
W, and to calculate its shape at the MIT. Our results suggest that the function Q n {e) shows 
critical behavior near W c . In Fig. [| we show Q n (e) corresponding to the critical point. 



Calculations for different system sizes at W c yield almost the same set of distributions Fig. 2 
independently of the number of levels n = 0,1,2,.... The width of Q n (e) at the MIT is 
larger than in the metallic regime, since the level repulsion becomes less important. We 
note that for small fluctuations (8n <C s) hiQ n (e) is quadratic with respect to both e and 
n, similarly to ([!]). 

Let us consider the quantity /„ = J °° Q n (e)de. For the Gaussian orthogonal ensemble 
it is known from previous numerical simulations || that I w 0.643, I\ ~ 0.922, and I n 
converges to unity in the limit n — > oo. For the Poissonian process (0) it is easy to see that 
/„ = 1 for all n. At the MIT, when W = 16.5, we obtained the following values of the 
integral: J C w 0.714, Jf w 0.960, J 2 C w 0.994, J 3 C w 0.998, and also lim^^ I c n = 1. The latter 
set of constants is L-independent and characterizes the critical level statistics. We emphasize 
that I n changes with W around the critical disorder W c according to a one-parameter scaling 
law, I n (W, L) = + f(L/£(W)), where the function f(x) can be linearized in the vicinity 
of the critical point, f(L,W) oc (W — W c )L 1 ' h '. This allows determination of the critical 
exponent v of the correlation length. Note that in searching for v it is not necessary to 
choose a certain spacing s Q as was done earlier when using P(s). [|IT1 ,|T2] 



IV. VARIANCE OF THE LEVEL NUMBER 

The width of the distribution Q n (e) describes the rigidity of the spectrum, which is 
defined by the variance of the level number, ([5n(e)} 2 ) = S.^Lo( n — ( n )) 2 Qn(£), where (n) = e 
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is the average number of levels in a given interval e. One can ask how dependence of the 
variance ([5n(e)} 2 ) on the average (11(e)) changes when the delocalized states transform into 
the localized ones. 

It follows from the random-matrix theory that in the metallic regime the variance 
is given by the Dyson formula, 



(S 2 n) M = \(\n(n) + C), 1 « (n) < hD/L 2 , (4) 



where C ~ 2.18 and D is the diffusion constant. By increasing the disorder W the fluc- 
tuations of 5n are increased due to a weakening of level correlations. In the strongly lo- 
calized regime the levels are not correlated, hence (S 2 n)i = (n). This is much larger than 
(8 2 n) M . Exactly at the MIT, the variance depends linearly on the average level number, fl5| 
(5 2 n) c = tt(n), as in the insulating limit. However the factor k is less than unity. 

We calculated the dependence of the ratio (5n 2 ) / (n) on the average number of levels 
within a given interval e for different lattice sizes L at W c = 16.5 as shown in Fig. 13. 



Obviously the ratio is not sensitive to a change of L. Such behavior results from the above- Fig. 3 
mentioned universality of the probability density Q n (e). One observes that the amplitude 
of the relative fluctuations of the number of levels (5 2 n) / (n) decreases very weakly with 
the energy, approaching the constant limit k ~ 0.27. Thus, our results are consistent with 
the suggestion of proportionality between the variance and the average level number at 
large (n), [15| but deviate from the recently proposed power law, (5 2 n) / (n) oc (n) _1 / 3v . 



However, this power law can be considered in our calculations as the next higher-order 
correction to the leading constant term. 

It was earlier conjectured [|15|] that Q n {^) at the MIT is normally distributed oc 



exp[— (e — n) 2 /2(5 2 n) c ] for fixed e 3> n. After substituting ke for the variance and dou- 
bly differentiating, for n = one recovers the Poissonian asymptotic form of the critical 
P(e) oc exp(— e/2k). 
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V. CONCLUSIONS 



We have discussed the statistics of the spectral fluctuations in the 3D systems at the 
disorder-induced metal-insulator transition. By directly diagonalizing the Anderson Hamil- 
tonian the nearest-neighbor level spacing distribution P(s) and the probability Q n (e) of 
finding exactly n eigenvalues in an energy interval of width e were numerically calculated. 
We confirmed that the critical P(s) is L-independent and has a Poissonian like asymptotic 
form at large s, but with stronger decay than that in the insulating regime. One expects 
that the critical Q n (e) which describes the statistical fluctuations in the discrete spectrum 
exhibits also complete scale-invariance in analogy with the universality of the critical P(s). 
Our results show that the distributions Q n (e) at the MIT are well described by Gaussians 
both for small fluctuations Sn <C e and also when n <C e, as in the metallic limit. We also 
analyzed the variance of the number of states (S 2 n) in a given interval as a function of mean 
level number (n). At the critical point this function was found to be independent of the 
system size and to have a leading linear term, (5 2 n) c = 0.27(n) for n >> 1. To summarize, 
our results indicate that the statistics of the energy levels near the Anderson transition con- 
stitute an infinite set of scaling variables I n which can be used to characterize the mobility 
edge by one-parameter scaling laws. 
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FIGURES 

FIG. 1. Level spacing distribution P(s) at the MIT for various system sizes L. Curves: Wigner 
and Poisson distributions. Inset shows the tail of P(s). Straight line is best fit. Dotted line: the 
power-law asymptotic || lnP(s) oc s l+l ^ u . 

FIG. 2. Critical probability distribution Q n {e) for different n. 

FIG. 3. Variance of the level number (5n 2 ) as a function of the average (n) for various L. 
Solid line: Dyson result (Eh. 
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